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Abstract. We will show how to get Burkholder's function from |Bul| by using 
Monge- Ampere equation. This method is quite different from those in the series 
of Burkholder's papers [Bul| - [Bu7| . 

1. Introduction 

Bellman function method in Harmonic Analysis was introduced by Burkholder 
for finding the norm in IP of the Martingale transform. Later it became clear that 
the scope of the method is quite wide. 

The technique, originated in Burkholder's papers [Bulj -[Bu7j. can be credited 
for helping to solve several old Harmonic Analysis problems and for unifying ap- 
proach to many others. In the first category one would name the (sharp weighted) 
estimates of such classical operators as the Ahlfors-Beurling transform (Banuelos- 
Wang [BaWal], Banuelos-Janakiraman [BaJalj, Banuelos-Mendez [BaMH], Na- 
zarov-Volberg |NVlj . Petermichl-Volberg |PVj . Dragicevic-Volberg |DV2j ) and 
the Hilbert and Riesz transforms (Petermichl |Plj . |P2j ). In the second category 
one can name all kind of dimension free estimates of weighted and unweighted 
Riesz transforms (see a vast literature in [DV1] - [DV3| ). Roughly, Bellman func- 
tion method makes apparent the hidden scaling properties of a given Harmonic 
Analysis problem. Conversely, given a Harmonic Analysis problem with certain 
scaling properties one can (formally) associate with is a non-linear PDE, the so- 
called Bellman equation of the problem. 

Let us recall to the reader that in the series of papers |Bu| . |Bul] - [Bu7] Donald 
Burkholder investigated Martingale transform and gave the sharp bounds on this 
operator in various settings-but by similar methods. The methods were so novel 
and powerful that the influence of these articles will be felt for many years to come. 
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The novelty was a key. One of the leading mathematician working in the domain of 
Harmonic Analysis told the second author that these papers of Burkholder "spin 
his head". In the book of Daniel Strook [Str] many pages are devoted to the 
technique developed by Burkholder in the abovementioned series of papers, and 
the reader can sense the same feeling. It is explained in [Str] that the simplest way 
to understand the sharp estimates of Martingale transform obtained by Burkholder 
is to operate with one of the so-called Burkholder's function: 

u p (x,y) = p(l - ^) p -\\y\ - (p* - l)\x\)(\x\ + lyir 1 , (1) 

here p* := max(p, r^j), 1 < p < oo . 

However, the main question is of course how to get this function? Where did 
it come from? These questions are asked in [Str] as well. Of course, Burkholder 
explains in many details the way this function (and several of its relatives) are 
obtained. It is almost (but not quite) the least bi-concave majorant of function 

\y\P~l - (p* _ l)P|a;|P . (2) 

It is obtained by solving a certain PDE and performing certain manipulations with 
the solution after that. The reader will find much more about u p after reading this 
article, in particular in Section [6J 

But it seems like the same questions persist even after this explanation. And 
a new question can appear: how wide is the applicability of the technique that 
Burkholder elaborated in [Bulj - [Bu7j ? There is a vague feeling that the area of 
applicability is quite wide. To make this feeling more precise one should look at 
the function above closer and see that it is a creature from another universe, which, 
initially, does not have too much in common with Harmonic Analysis. Burkholder 
function is a natural dweller of the area called Stochastic Optimal Control. It 
is a solution of a corresponding Bellman equation (or a dynamic programming 
equation) but in the setting, when the differential equations subject to control are 
not the usual ones. They are stochastic differential equations. The reader can 
find some notes on this in [VoEcoIej . |JNTV2j . |VaVo| . |VaVo2j . [SlSt] , These notes 
explain why Stochastic Optimal Control is the right tool to work with a certain class 
of Harmonic Analysis problems. On the other hand, Stochastic Optimal Control 
problems generically can be reduced to solving a so-called Bellman PDE (and 
proving the so-called "verification theorems", but this is a second task). Bellman 



BURKHOLDER'S FUNCTION 



3 



PDEs belong to the class of fully non-linear PDEs. Often they are PDEs of Monge- 
Ampere type. In the present article we would like to show the reader how to obtain 
Burkhloder functions (the one above and others from [Bul| - [Bu7| ) by reducing the 
search for them to solving certain Monge-Ampere equations. The scope of the 
application of the methods of Stochastic Optimal Control to Harmonic Analysis 
proved to be quite large. After Burkholder the first systematic application of this 
technique appeared in 1995 in the first preprint version of [NTVlj . It was vastly 
developed in [NT] and in (now) numerous papers that followed. A small part of 
this literature can be found in the bibliography below. 



2. Notations and definitions 

We shall say that an interval / and a pair of positive numbers af such that 
a + + a~ = 1 generate a pair of subintervals I + and I~ if {1^1 = ctf\I\ means 
the length of I) and I = I~ U / + . For a given interval J the symbol J = J {a) will 
denote the families of subintervals of J such that 

• J G J] 

• if J G J then J± G J . 

For a special choice if all of = | we get the dyadic family J = T>. Every family 
J has its own set of Haar functions: 



VI G J h x (t) = < 



if t el- 
ate u 



If the family J is such that that the maximal length of the interval of n-th gener- 
ation (i.e., after splitting the initial interval J into 2 n parts) tends to as n — > oo, 
the Haar family forms an orthonormal basis in the space L 2 (J) Q {const}. 



For a function / G L l (I) the symbol means the average of / over the interval 



I: 



if) j 



f(t)dt . 



Definition. Fix a real p, 1 < p < oo, and let p' = p* = m&x{p,p'}. Introduce 
the following domain in M 3 : 



Cl = tt(p) = {x = (xi,x 2 ,x 3 ) : x 3 > 0, |xi| p <x 3 }. 
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For a fixed partition J of an interval J we define two function on this domain 
B max (a;) = B max (x;p) = sup{(|#| p ) A, 

f,9 

B min (x) = B min (x;p) = hrf{{\g\ p )j}, 

f,9 

where the supremum is taken over all functions /, g from L P (J) such that (/) 7 = xi, 
(g) j = X2, {\f\ p )j = %3, and \(f,hj)\ = \(g,hj)\. We shall refer to any such pair of 
functions /, g as to an admissible pair. When \(f,hj)\ = \ {g,hj)\ happens for all 
dyadic intervals inside J we call g a Martingale transform of /. We shall call 
B max (x) (and B m i n (x)) the Bellman functions of the problem of finding the best 
constant for the Martingale transform inequality: 

\(9)A<\(f)j\^ (\9\ P h <C(p)(\f\ P )j- (3) 

This best constant was found by Burkholder: C(p) = (p* — l) p , p* := max(p, ^j). 
Remark 1. It is amazing that there is no proof that would find this C(p) without 
finding the function of 3 variables B max (x) or some of its relatives (like, for example, 
u p from (H])). 

Remark 2. Burkholder proved that the functions B do not depend on the initial 
interval J and on a specific choice of its partition. Below we work only with dyadic 
partitions. 

Remark 3. In the case p = 2 the Bellman function are evident: 

B max ( 2; ) = B m i n (a;) = x\ + X3 — x\ . 



Indeed, since 



we have 



i = |J|x 3 = |J|x? + ^|(/, h x ] 



^\\9\\l = 4 + ^E\^h I )\ 2 

x 2 + nr 12 K/' = x 2 + ^3 - x\ 
2 



Define the following function on Ri. = {z = {zi,z 2 )'- zi > 0}: 

„, v \[zi-(jp*-iyz p i if^<( P *-i)^ 2 , 

F p (z 1 ,z 2 ) = < (4) 

( P {l-±) p - 1 (z 1 + Z 2 ) p - 1 [z 1 -( P *-l)z 2 \, if Zl > ( P *-1)Z 2 . 
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Note for for p = 2 the expressions above are reduced to F 2 (zi, z 2 ) = Z\ — z\- 



3. The main result 
Now we are ready to state the main result: 

1. 1 

Theorem 1. The equation -F p (|a:i|, |x2|) = -Pp( x f,B p ) determines implicitly the 
function B = T5 m [ n (x;p) and the equation F p (\x 2 \,\xi\) = F P (B P ,x%) determines 
implicitly the function B = B max (x;p). 

Remark. The reader can take a look at formulae (5.23)— (5.27) on page 660 of |Bul] 
and recognize that this is how Burkholder describes B max . The same is true for 

B m in- 

4. How to find Bellman functions 

We start from deducing the main inequality for Bellman functions. Introduce 
new variables y\ = \{x 2 + si), 2/2 = \{x 2 — and 2/3 = x 3 . In terms of the new 
variables we define a function M, 

M(yi, 1/2,2/3) = B(xi,x 2 ,x 3 ) = B(yi - y 2 ,yi + 2/2,2/3) , 
on the domain 

2 = {y = (2/1,2/2,2/3) : 2/3 > 0, \yi - y 2 \ p < 2/3} • 

Since the point of the boundary x% = \xi\ p (2/3 = \y± — 2/2| p ) occurs for the only 
constant test function f = x\ (and therefore then g = x 2 is a constant function as 
well), we have 

B(xt,x 2 , \xi\ p ) = \x 2 \ p , 

or 

M(2/i,2/2, 12/1 -2/2| p ) = I2/1 +2/2| p - (5) 
Note that the function B is even with respect of x\ and x 2 , i.e., 

B(xi,x 2 ,x 3 ) = B(-xi,x 2 ,x 3 ) = B(xi, -x 2 ,x 3 ) . 

It follows from the definition of B if we consider the test functions / = — / for the 
first equality and g = —g for the second one. For the function M this means that 
we have the symmetry with respect to the lines 2/1 = ±2/2 

M(2/i,2/2,2/3) = M(2/2,2/i,2/3) = M(-yi, -y 2 , 3/3) • (6) 
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Therefore, it is sufficient to find the function B in the domain 

ft + = fi + (p) = {x = (xi,x 2 ,x 3 ) : Xi > 0, |xi| p <x 3 }, (7) 
or the function M in the domain 

H+ = {y = (yi, 2/2, 2/3) : yi > 0, -yi < 2/2 < yi, (yi - y2) p < ^3} • (8) 

Then we get the solution in the whole domain by putting 

B(x 1 ,x 2 ,x 3 ) = B(| xi\, I x 2 \,x 3 ) . 

Due to the symmetry (|6j) we have the following boundary conditions on the "new 
part" of the boundary dE + : 

dM dM 

- — = - — on the hyperplane y 2 = yi , 
oyi dy 2 

dM dM , , U 

— — = —-7: — on the hyperplane y 2 = —y± . 
oyi oy 2 

If we consider the family of test functions / = rf, g = rg together with / and 
g we come to the following homogeneity condition 

B(tx 1 ,tx 2 ,t p x 3 ) = t p B(x 1 ,x 2 ,x 3 ) , 

or 

M(ry 1 ,Ty 2 ,T p y 3 ) = T p M(y 1 ,y 2 ,y 3 ) . 

We shall use this property in the following form: take derivative with respect to r 
and put r = 1 

dM dM dM 

yi-s Ky2^ \-py3-K— =pM(y u y 2 ,y3) ■ (10) 

oyi dy 2 dy 3 

Let us fix two points x^ G $7 such that \xf — x^[\ = \x^ — x^l, for the cor- 
responding points y^ G 5 this means that either yj*~ = y7/ ; or y^ = y 2 ■ Then 
for an arbitrarily small number e > by the definition of the Bellman function 
B = B max there exist two couples of test functions / and g^ on the intervals 
such that (/% = xf, (y±) j± = xf , fl/^)^ = xf, and fly^)^ > B(x±) - e. 
On the interval I = I + U I~ we define a pair of test functions / and g as follows 
f\I^ = / ± , y|/ ± = y 1 * 1 . This is a pair of test functions that corresponds to the point 
x = a + x + +a~x~ , where a ± = |.Z' ± |/|.Z'|, because the property \x^— xj~| = Ix^—x^l 
means |(/, h T )\ = \(g,hj)\. This yields 

B(s) > (|y| p ), = a+(|y+| p )+ + a-(|y-| p )7 > a+B(x+) + arB(ar) - e . 
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Since e is arbitrary we conclude 

B(x) > a+B(x + ) + a~B(x~). (11) 
For the function B = B m i n we can get in a similar way 

B(s) < a + B(x + ) + a~B(x"~). (12) 

Recall that this is not quite concavity (convexity) condition, because we have 
the restriction \xf — x~^\ = \x^ — x^l- But in terms of the function M 

M max (y) > a+M max (y + ) + a~M m! , x (y-) , 
M min (y) < a + M min (y+) + cTM^y-) , 

when either yi = = Vi , or y2 = y£ = y^i we indeed have the concavity 
(convexity) of the function M with respect to yi-, 2/3 under a fixed yi, and with 
respect to yi, y% under a fixed y%- 

Since the domain is convex, under the assumption that the function B are suf- 
ficiently smooth these conditions of concavity (convexity) are equivalent to the 
differential inequalities 

(Wy iyi My iy3 \ ^ fMy 2 y 2 My^yA ^ ^ V?/ € H , (13) 

d 2 M 

for M = M max (here My iBj stand for the partial derivatives — — - — ) and 

dyiuyj 

fMy iyi M yi y 3 \ Q ^ fMy 2 y 2 My^ \ Q ^ Vtf € 3 , (14) 

for M = M min . 

Extremal properties of the Bellman function requires for one of matrices in f)13|) 
and (fT4"|) to be degenerated. So we arrive at the Monge-Ampere equation: 

M M ,M M3 = (M ym ) 2 (15) 

either for i = 1 or for i = 2. To find a candidate M for the role of the true Bellman 
function M we shall solve this equation. After finding this solution we shall prove 
that M = M. 

The method of solving homogeneous Monge-Ampere equation is described, for 
example, in |VaVo| . [VaVo2] . [SIStj . In particular we know that the solution of the 
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Monge-Ampere equation has to be of the form 

M = kVi + t 3 y 3 + t , (16) 

where if. = M yk , k = 1,2,3. The solution M is linear along the lines (let us call 
them extremal trajectories) 

yidU + y 3 dt 3 +dt o = 0. (17) 

One of the ends of the extremal trajectory has to be a point on the boundary 
2/3 = 1 2/1 — 2/2 | p , where constant functions are the only test functions corresponding 
to these points. Denote this point by U = (y%,u, (yi — u) p ). Note that we write 
(yi — u) p instead of |yi — u\ p because the domain H + is under consideration. For 
the second end of the extremal trajectory we have four possibilities 

1) it belongs to the same boundary y 3 = (yi — y2) p ; 

2) it is at infinity (yi,y2,+oo), i.e., the extremal lines goes parallel to the 
y 3 -axis; 

3) it belongs to the boundary yi = y\\ 

4) it belongs to the boundary y2 = —y\. 

The first possibility gives us no solution. Namely, we have the following 

Proposition. If p ^ 2, then the function B max cannot be equal to B{x) = M{y), 
where M is the solution of the Monge-Ampere equation (|15p such that one of its 
extremal trajectory is of type 1) above. The same claim holds for B m ; n . 

Proof. To check this it is sufficient to verify that the test functions of the type 
a + f3hj (t) cannot be an extremal function of our problem with the only exception 
of p = 2, when the situation is trivial: B max (x) = B m ; n (x) = x 3 + x^ — xf, and any 
pair of test function is extremal. We will show that the Bellman functions being 
solution of the homogeneous Monge-Ampere equation cannot be linear on a chord 
connecting two points x^ on the boundary <9$7, i.e. such a chord cannot 
be an extremal trajectory of our Monge-Ampere equation. 

We assume now that two points x^ S Q + such that \x"T — x^\ = \x^ — x^\ 
are on the boundary xf = (xf) p and x = ^(x + + x~). We need to show that 
\{{ x t) P + ( x 2 ) P ) can be the value neither of B max (a;) nor of B m i n (x). 

Without lost of generality we may assume that xf > x\. Let us denote a := 
^(xf — x±) then 

xf = xi db a xf = X2 db a a , 
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where a = ±1 depending on the direction of our chord: it can be either in the 
plane x\ — x% = const (and then a = 1) or in the plane x\ + x% = const (and then 
cr = —1). The pair of the test functions /, g on / = [0, 1] that gives the value 

Ma) =: <|#>, = \(\xt\ p + I * 2 T) = ^(1*2 + o|* + |*2 - 

is 

/ = xi + a/i 7 (7 = £2 + • 

First of all we assume that x^ 6 but xf ^ i.e. xi ±a > 0, X2±a > 0, because 
if one of xf is zero we are in the cases either 3) or 4) listed before Proposition. 
Our aim will be to find another pair of test functions f,g corresponding to the 
same point x, but with either bigger than A{a) (and then A(a) cannot be 

the value of B max (x)) or less than A{a) (and then A(a) cannot be the value of 
B min (x)). 

Let as make here two remarks. First, we see that the expression A{a) does not 
depend on the direction a. Therefore, the construction of the desired f,g ensures 
us that the point x cannot be the center of an extremal trajectory with two ends 
on dO, in any direction a = ±1. Secondly, we note that it is not obligatory to look 
for /, g for all a £ (0, max{xi, X2}), it sufficient to do this for small values of a/x2- 

Indeed, suppose that the chord L = [x~,x + ] represents an extremal trajectory 
of the corresponding Monge-Ampere equation. Let us consider the "crescent" 
between the chord L and the boundary d£l. It should be filled in by chords on which 
B(x) are linear (this is the property of the solutions of the homogeneous Monge- 
Ampere equation expressed in Pogorelov's theorem, see [Fog]). Among these chords 
we can take one (say L = [x~,x + ]) of arbitrarily small length (arbitrarily small 
value of a/x2)- Therefore, we can work with a new point x and new chord L: would 
we show that L cannot be an extremal trajectory, the chord L could not be one 
either. 

To construct the desired pair /, g we use the following family of functions 
equal to 1 on [0, 1/2 — s] U [1 — s, 1] and to —1 on (1/2 — s, 1 — s). Let us note that 
all (fi s have the same distribution function as hj, and 

(<t>s,hj) = 1 - 4s. 

We are interested in (f> := (fti/g. Then 
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and ip = 4> — hj has 

(il>,hj) = --, (tp,hj) = (<t>,hj) 

for all other dyadic J. So tp is a martingale transform of 4> (it is equal to on 
[0,3/8] U [1/2,7/8] and to ±2 on two intervals (7/8,1), (3/8,1/2)) and we can 
examine the pair 

/ = Xi + acf), g = x 2 + aip . 

Since / and / have the same distribution function, we have ( \ f\ p ) I = { |/| p ) 7 = £3, 
i.e., f,g is a pair of test functions corresponding to the same point x. 
To investigate the difference ( \g\ p ) I — ( \g\ p )j we use the function 

Ap(a) := ^((1 + 2af + (1 - 2a)*) + \ - ((1 + a) p + (1 - af) . 

Since 

Ap(a) = ^ P (p - l)(p - 2)(p - 3)a 4 + 0(a 6 ) , 

we have A p (a) > for small aifl<p<2orp>3 and X p (a) < for small a if 
2 < p < 3. Recall that 

z 2 < t < § 
x 2 - 2a § < i < \ 



x 2 \<t<l 
x 2 + 2a l<t<l, 



therefore, ( \g\ p ) I = |((x 2 + 2a) p + (x 2 - 2a) p ) + fx^ and 

<l#),-(l#>,=^A(f ). 

x 2 

For small a we have a desired example for B max if p G (1,2) U (3, 00) (because 
A p > and for B m i n if p € (2, 3) (because A p < 0). 

Now we interchange in a sense the roles of / and g: instead of a we take a new 
parameter, say a, and put 

/ = x\ + dip, g = X2 + a<fi. 

Since we have now 



(i/a-(i/a = ^A(f ) 



and the function 



t h+ ( |x + tVr)/ = -(k - 2t| p + |x + 2i| p ) + ^|x| p 
8 4 
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is increasing in t > from \x\ p till infinity, we can find a, a > a for p £ (2, 3) and 
a < a for p G (1, 2) U (3, oo), such that 

(m p >, = <l/r°>,=z3. 

For this a we get the desired pair of test function, because the function 

(\x + t4>\P) I = ^(\x-t\ p + \x + t\ p ) 

is also increasing in t > 0, and therefore, we have 

( \~g\n, = (\x 2 + ~a^) 1 > ( |* 2 + ~acj>\ p ) 1 = ( 1*2 + ^H, = < M P >i 

if p € (2, 3) and the opposite inequality if p € (1, 2) U (3, oo) (and a < a). 

This construction failed for p = 3 because A3 (a) = for all a 6 (0, ^). To avoid 
this difficulty we modify the function ip, namely, we take ip = <p + h i+ /y/2, i.e., 



1, 


o<;<! 


-1, 


f <*< ^ 


< o, 


i<t<| 


-2, 


i<<<! 




|<t<i 



The function ip is a martingale transform of </>, since 
Now we put 

/ = x\ + a0, <7 = *2 i . 
As before we have ( |/| 3 ) 7 = \ {{x - af + {x + a) 3 ) = ( |/| 3 ) 7 , but 

( l^l 3 ), = x% + 6x 2 a 2 T^a 3 = ( \g\ 3 ) I T ^a 3 , 

therefore, by choosing the sign in the definition of g we are able to increase as well 
to decrease the value ( |g| 3 ) 7 , hence it is neither the value of B max (x) nor B m ; n (x). 
Proposition is completely proved. □ 

Now we check the second possibility among the possibilities l)-4) listed right be- 
fore the Proposition. Since the extremal line is parallel to the 2/3-axis, the Bellman 
function has to be of the form 



M (y) = A {yi, V2) + C(y 1 ,y 2 )y3 ■ 
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Any pair of inequalities both ([TBI and (fT4"l) implies M ym My 3 y 3 — {M ym ) 2 > 0. 
Since M y , iyi = 0, this yields M ym = = 0, i.e., C is a constant. From the 
boundary condition ([5]) we get 

A(yi,y 2 ) + C(yi - y 2 ) p = (yi + y 2 ) p , 



whence 



and 



or 



A(yi , 2/2 ) = (yi + y 2 ) p - C(yi - y 2 ) p , 
M(y) = (yi + y 2 f + C(y 3 - ( Vl - y 2 f) , (18) 



B(x) = x p 2 + C(x 3 - x{) . (19) 

Let us note that this solution cannot satisfy necessary conditions in the whole 

domain H + except the case p = 2. The constant C must be positive (otherwise 

the extremal lines cannot tend to infinity along y3-axes, because M must be a 

nonnegative function). Therefore, the straight line 

i i 
yi + y2 = Cp- 2 (yi -y 2 ), orx 2 = 0- 2 yi 

splits S + in two subdomains, in one of which the derivatives 

d 2 M d 2 M . / 2 2 \ 

~dyf = ~dyY = ~ ^ " " ^ } ) 

is positive (i.e., it could be a candidate for B m j n ), and in another one is negative 

(i.e., it could be a candidate for B max ). 

Thus, this simple solution cannot give us the whole Bellman function and we 
need to continue the consideration of the possibilities 3) and 4) (listed right be- 
fore the Proposition. Till now we have not fixed which of two matrices in (|13p or 
in (fl~4"|) is degenerated, i.e., what is i in the Monge-Ampere equation (fl~5j) . because 
for the vertical extremal lines both these equations are fulfilled. Now, when con- 
sidering possibility 3) or 4), we need to investigate separately both Monge-Ampere 
equations (fT5|) . We shall refer to these cases as 3j) and 4j). 

Let us start with simultaneous consideration of the cases 3i) and 4i) (we recall 
that this means that y 2 is fixed). We look for a function 



M = tm + t 3 y 3 + t 
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on the domain H+, which is linear along the extremal lines 

yidh + y 3 dt 3 + dt = . 

Now one end point of our extremal line V = (v,y 2 , (v — 2/2 ) p ) belongs to the 
boundary y 3 = \yi — y 2 \ p and the second end point W = (|2/2|j 2/2, w) is on the 
boundary y\ = \y 2 \-, where we have boundary condition ([9]). Due to the symme- 
try ©, on the boundary y% = y2 (this means that our fixed y 2 > 0) we have 

dM dM 



and 



dy2 dyi 
dM dM 



dy 2 dyi 

on the boundary y\ = —y 2 (this means that our fixed y 2 < 0). In both cases 

dM dM 
dy 2 dyi 

and therefore (fTUj) and (fTBT) imply 

2*i 1 3/2 1 +pwt 3 = pM(W) = ph\y 2 \ + pwt 3 + pt , 

whence 

,2 , 

*0 = ( 1)*1 2/2 • 

This gives the formula for *q(*i) (remember that y 2 is fixed as we consider the cases 
3i), 4i) now). Thus, we get 



M(y) 



yi + (- - l)|y 2 | 



*1 + 2/3*3 • 



(20) 



Since d*o = (f — 1 ) 1 2/2 1 <**i> the equation of the extremal trajectories (fl~7l) takes 
the form 



2/1 + (- - III2/2I 



and we can rewrite (1201) as follows 



dti + y 3 dt 3 = 0, 



(21) 



We see that the expression M(y)/y 3 is constant along the trajectory and we can 
find it evaluating at the point V, where the boundary condition (|5|) is known: 



M(y) 



y + y2 

V - 2/2 



2/3, 



(22) 
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where v = f(2/i) 2/2> 2/3) satisfies the following equation: 

1/9 I 

(23) 



2/x -h (| — 1) I2/2I _ u + 



2/3 - 2/2) P 

because the point V = (v,t/2,(v — IJ2Y) is on the extremal line (|2~T1) . We even 
shall not check under what conditions equation (I23p has a solution and when it is 
unique. Later we show that in any case the function M we have found cannot be 
the Bellman function we are interested in, because neither condition (j!3p nor f)14[) 
can be fulfilled: the matrix {M yiy . }ij=2,3 is neither negative definite nor positive 
definite. We postpone this verification, because the calculation of the sign of the 
Hessian matrices is the same for this solution and another solution of the Monge- 
Ampere equation that supplies us with the true Bellman function. And these 
calculations will be made simultaneously a bit later. And now we only rewrite our 
solution in an implicit form more convenient for calculation. 
We introduce 

then @ yields 

v = ^V2. (25) 
uj — 1 

Since v > (in fact, recall that we consider now only y: 2/1 > I2/2I domain now, and 
that v is just the first coordinate of the point V = (v , 2/2, (v — 2/2 ) p in this domain), 
we have 

sign 2/2 = sign(ci; — 1) . (26) 
After substitution of ([25]) in ([23]) we get 



22/2 
uj - 1 



.2 



2/1 + ( l)|2/2| =2/3 72/2 + ( 1) I2/2 

p - La; — 1 p 



w + 1 ,2 



V 
or 

2 p |2/ 2 | p - 1 [rai + (2 - p)\V2\] = 2/3^ " 1| P_1 [(w + 1)P + (2 - p)|w - 1|] 

For the case 3i) we have 2/2 > (i.e., £2 > ^lj we look for uj > 1 or B > 2/3) and 
the latter equation can be rewritten in the initial coordinates as follows 

{x 2 - xif- 1 [x 2 + (p - l)xi] = (B* - xlf- 1 [B* +{p- 1)4] • 

For the case 4i) we have 2/2 < (i.e., X2 < xi, we look for uj < 1 or B < 2/3) and 
the equation takes the form 

On - x 2 ) p ' 1 [xi + {p- l)x 2 ] = (xf - B^f- 1 [x\ + (p- 1)B'] 
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Introduce the following function 

G( Zl ,z 2 ) = (z x + z 2 f- 1 [zi - (p - l)z 2 ] 
denned on the half-plane Z\ + z 2 > 0. Then in the case 3i) we have the relation 

G{x 2 ,-x x ) = G{B\-x\), (27) 

or 

G(y 2 + 2/1,2/2 - yi) = y 3 G(uj,-l) . 

In the case 4i) we have 

G(x 1 ,-x 2 ) = G(x*,-B*), (28) 



or 



G(yi - y 2 , -yi - 2/2) = 2/3^(1, -u) 



Now we have to consider the Monge-Ampere equation (I15p in the cases 32) 
and A 2 ). This means that we fix y\ now. Let us begin with the cases 32), when an 
extremal line starts at a point U = (yi,u, (2/1 — u) p ) on our parabola and ends at 
a point W = (yi,yi,w). Again, the symmetry condition at the point W is 

dM _dM _ 
dyi dy 2 

and the homogeneity condition (|10p plus condition (|16p at W yield 
2yit 2 + pwt 3 = pM(W) = pyit 2 + pwt 3 + pt , 

whence 



*o = ( 1)2/1*2 , 

p ' 



and therefore 



M(y) 



2/2 + (- - IJ2/1 



*2 + 2/3*3 • 



(29) 



Since dto = (| — l^jyidt 2 , the equation of the extremal trajectories takes the 
form 



2/2 + (- - 1)2/1 



and we can rewrite (1291) as follows 



d*2 + 2/3 rf*3 = , 



(30) 



M(2/) = [t 3 -t 2 —^)y 3 . 
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Again, from here the expression M{y)/y^ is constant along the trajectory and we 
can find it evaluating at the point U, where we the boundary condition ([5]) is 
known: 

M(y)=( y -^Xys, (31) 
where u = u(j/i, 2/2, 2/3) can be found from ([30]) : 

2/2+ (| -1)2/1 = ^+(|- 1 )2/i 

2/3 (2/1 - u) p 

We see that if our extremal line starts at point U = (yi,u, (y\ — u) p on our 
parabola u = — (| — l)j/i, then y2 = — (| — 1)2/1 = u = const, i.e., it is a line 
parallel to the X3-axes. This means that no extremal line that ends at the points of 
the boundary y\ = y2 can intersect the plane yi = — (| — l)yi. This follows from 
the property that extremal trajectories do not intersect. Therefore, the starting 




2/2 = yi 



yi 

V2 = -(I - l)yi 



V2 = -J/i 



Figure 1. Acceptable sector for the case 32). 



points U with u < — (| — l) y\ cannot be acceptable for the case under consideration 
(since these trajectories do not intersect the plane yi = — (| — they cannot 
have the second end point on 2/2 = 2/1, see Fig. [T|). 

Let us check that equation ([32]) has exactly one solution u = u(yi, 2/2, 2/3) i n the 
sector — (- — l)yi < 2/2 < 2/1- Indeed, the function 



« ^ 2/3 



n + ( ljyi 



(2/1 



2/2 + (- - 1)2/1 
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is monotonously increasing for u < y\ and it has the negative value —(-niV' 



2/2 + 



1)2/1 



at the point u 



l)yi and the positive value |yi2/3 at the point 



u = yi. 

Now we rewrite the solution (I3ip in an implicit form using notations ([2 



M {y)\ v 

J/3 



From (1311) we have 



u 



CO - 1 
w + 1 



2/1 , 



(33) 



therefore, from (|32p we obtain 

2^(0; + l)^ 1 [p(w - 1) + (2 - p)(u + 1)] = y?- 1 [py 2 + (2 - p) yi ] 

or 

2-P+i 1 fe(o> + l)*" 1 ^ - p + 1) = y?- 1 [py 2 + (2 - p) yi ] , 
which is (using again notations ([24")) : w := P ) 

(5* + xly- 1 [B* - (p - = (xi + x 2 ) p_1 [x 2 - (p - l)xi] . 

In terms of function G this can be rewritten as follows 

G(x 2 ,x 1 ) = G(B*,xj), 



or 



G(yi + 2/2, 2/i - 2/2) = 2/3G(w, 1) 



-2yit 2 + pwt 3 = pM(W) = -py\t 2 + pwt 3 + pt , 



It remains to examine the possibility A 2 ). Assume that an extremal line starts 
at a point U = (yi,u, {y\ — u) p ) and ends at a point W = (y%, —y%,w). Again, the 
homogeneity property (fTUj) at the point W and the symmetry = — = —t 2 
yield 

whence 
and therefore 

Since dto = 
form 



1- - 2/1*2, 



M{y)= y 2 + (l--) yi 
L p 



*2 + 2/3*3 • 



(34) 



(1 — ^)yidt 2 , the equation of the extremal trajectories takes the 



2/2 + (1 ) 2/1 d*2 + 2/3 dt 3 = , 

p J 



(35) 
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and we can rewrite (1341) as follows 

M{y) 



, , dt 3 \ 



Again, the expression M(y) /y% is constant along the trajectory and from the bound- 
ary condition ([SD we get the same expression 



M{y) 



yi + u 



2/3 



(36) 



.yi - u , 

Now u = u(yx, 2/2, 2/3) is a solution of the equation 

m - (| - _ (| - i)&fi 
y3 (2/1 - 

that we get from ([35]) . As before, we get trajectories ending at the plane y 2 = — y\ 
not in the whole domain H+, but only in the sector — y\ < y 2 < (| — l)yi (see Fig. [2]), 



(37) 



and equation (|37p has a unique solution for every point from this sector. As before, 



J/2 = 2/1 




Figure 2. Acceptable sector for the case 42). 



relation (J33]) allows us to rewrite the equation of extremal trajectories (137p as an 
implicite expression for ui (and hence for M): 

G(x 1 ,x 2 ) = G(xg,B'), 

or 

G{y\ - ?/2, 2/1 + JV2) = 273^(1, u) . 
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Now we start the verification which of the obtained solutions satisfies condi- 
tions (fT3l) or (TB|) . We need to calculate Di := My iyi M y3 y 3 — My. y3 , i = 1, 2, in four 

cases 

3i) G{yi + y 2 , -yi + y 2 ) = y?,G{uj, -1); 
4i) G{yi - y 2 , -yi - y 2 ) = 2/3^(1, -u); 

3 2 ) G(y 1 +y 2 ,y 1 -y 2 ) = y 3 G(u;,l); (38) 

h)G(y 1 -y 2 ,y 1 + y 2 ) = y 3 G(l,u;), (39) 
where M = y^oj p . In all situations we have a relation of the form 

2/3 

Till some moment in the future we will not specify the expression for $ and H, 
as well as for their derivatives, and plug in the specific expression only in the final 
result after numerous cancellation. In particular, we introduce 

Rl=Rl ^y-L and R 2= R 2 ( y: ).- R > 1 = - i ?L, 

We would like to mention here that this idea, allowing us to make calculation 
shorter, is taken from the original paper of Burkholder |Bul| . 
First of all we calculate the partial derivatives of lo: 

, H R\H 

^=-12 "y 3 =— -r> 

2/3 S/3 

§ uj v = — 1 =>• uj v = - = , 1 = 1,2. 

2/3 ' 2/3 2/3 

Here and further we shall use notation H' for any partial derivative H y ., i = 1,2. 
This cannot cause misunderstanding because only one i participate in calculation 
of Hessian determinants D{. Moreover, we shall not mention anymore that the 
index i can take two values either i = 1 or i = 2. 

vl vi vi 

R 2 cj y .H R X H' RxH' 
u y y — j — — ^—{R 2 H + y 3 J , 

yi yt yi 

2/3 2/3 2/3 



20 VASILY VASYUNIN AND ALEXANDER VOLBERG 

Now we pass to the calculation of derivatives of M = y^u p : 

M y 3 = PVi^ l ^y 3 + U p , 
M y . = py- i uj p ' 1 ujy i ; 

M y 3 y 3 = Py^ P ~ lu) y»yz + ^ p ' X ^y z + p(p- l)y 3 uj p ~ 2 uj 2 3 

= PUjP ~ 2 f lH2 [u;R 2 + {p- l)Ri] , (40) 

m 

M y 3 y, = PVzu^Uy^y. +puj p ~ 1 uj y% +p(p- l)y 3 u p ~ 2 Uy 3 Uy. 

pu p ~ 2 R 1 HH' . „ 

= t \°° R i + (P - l ) R i\ > 

2/3 



([^ J R 2 + (p-l) J R 1 ]( J H'0 2 +^3^' / ) • 



This yields 



2 2p-3n2 rr2 rr// 

A = My^yMy^ ~ M^ = ^ [ W tf 2 + - 1 ) ^ ] . (41) 

Notice, that ii 7 disappeared completely. 

Now we need to calculate second derivatives of 

H(yi,y 2 ) = G{aiyi + a 2 y 2 ,Piyi + #22/2), 
where aj,/3j = ±1. And 



dyf 



^" = 7TT2 ^"lyi + "22/2, /3iyi + /?2y2) 



= Gz A + Gz 2 z 2 ± 2G Zl 2 2 , 

where the "+" sign has to be taken if the coefficients in front of yi are equal and 
the "— " sign in the opposite case. 
The derivatives of G are simple: 

G Z1 = p( Zl + z 2 ) p ~ 2 [ Zl - (p - 2)z 2 ] , 
G Z2 = -p{p- l)z 2 ( Zl + z 2 ) p ~ 2 ; 
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G ZlZ2 = P (p - 1)(Z! + z 2 )P~ 3 [ Zl - (p - 3)z 2 ] , 
G Zl z 2 = -pip - l)(p - 2)z 2 (zi + z 2 )^ 3 , 
G % = -P(P " l)(*i + z 2 ) p ' 3 [zi + (p - l)z 2 ] • 

Note that (3^+6^= 2(7^ , and therefore, i?" = 4(7^ if oii = fa and H" = 
if Qj = — fa. The first case occurs for Hy 2 y 2 in cases 3i), 4i) and for i^ l2/l in cases 32), 
4 2 ). The second case occurs for Hy iVi in cases 3i), 4i) and for Hy 2 y 2 in cases 3 2 ), 4 2 ). 
In fact, we know that the equality D{ = has to be fulfilled in the cases 3j) and 4j), 
because it is just the Monge- Ampere equation we have been solving. 

So we have 



3i) 


Zl 


= VI +V2, 












Z2 


= -yi + 2/2, 


G ' 2 l 2 2 


= p{p- 


l)(p- 


2)(yi-2/ 2 )(22/ 2 r 3 , 


4i) 


Zl 


= yi - 2/2, 












Z2 


= -2/1 - 2/2, 


Gz l z 2 


= p{p- 


l)(p- 


2)(yi+2/2)(-2y 2 r 3 


3 2 ) 


Zl 


= 2/1+2/2, 












Z2 


= 2/1 - 2/2, 


Gz l Z 2 


= -p{p 


-l)(p 


-2)( yi -y 2 )(2yir 3 


4 2 ) 


Zl 


= 2/1 - 2/2, 












Z2 


= 2/1 +2/2, 


Gz l Z 2 


= -p{p 


-l)(p 


-2)(y 1 + y 2 )(2 2/1 ) p - 3 



In the first pair of cases we have signG zlZ2 = signi?" = sign(p — 2) and the 
opposite sign in the second pair of cases. In the first pair of cases we have that this 
sign is the sign of the Hessian determinant _D 2 up to the sign[u?i? 2 (p — l)Ri] (and 
D\ = identically); in the second pair of cases we have this sign is the sign of the 
Hessian determinant D 1 up to the sign[wi? 2 (p — l)R\] (and D 2 = identically). 

By the way, we call the attention of the reader to the fact, that, for example, in 
4i) above we have necessarily y 2 < (here y 2 is fixed and our extremal trajectories 
in the plane (2/1,2/3) here hit 2/1 = — 2/2 = 1 2/2 1 as we are always under restrictions 
—2/1 < 2/2 < 2/1, that is yi > \y 2 \), so (— 2y 2 ) p ~ 3 makes a perfect sense. The same 
type of observation holds for all other cases. 
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To complete the investigation of signDj we need to calculate the sign of the 
expression in the brackets in (|41l ): 



(42) 



ojR 2 + {p- l)Ri = R\\{p - l)& - 



(p 



4i) 



(p 



(P 



(P 



$(a; 
$'(w 

1)$' - £J$' 
$((J 

$'(a; 

1)$'-£J$' 
$((J 

$'(u; 
1)$' _ 



G(u;,-1), 

(w, -1) = P (cj - l) p - 2 (uj + P -2), 
G HH (u, -1) = p(p - l)( w - l) p " 3 (a; +P " 3), 

-Kp-l)(p-2)(w-l) p - 3 ; 

-G 22 (1, -w) = - P (p - l)w(l - w - l) p - 2 , 
G z ^ (1, -a,) = p(p - 1)(1 - ^) p " 3 [l - (p - 1)0;], 

-p(p-l)(p-2Ml-w)P- 3 ; 

G(u;,l), 

G^(w,l) =p(^ + l) p - 2 (^ -p + 2), 

G, A (w, 1) = p(p + l) p - 3 (a; - p + 3), 

-p(p-l)(p-2)(uj + l) p - 3 ; 

G(l,u), 

G Zi (l,u;) = -p(p-lMu J + iy- 2 , 

G ZxH (u>, -1) = -p(p - l)(w - l)f~ 3 [l + (p - l)w], 

-p(p-l)(p-2)a;(w + l) p - 3 . 



By the way, we call the attention of the reader to the fact, that, for example, in 
4i) above we have necessarily jte < 0, so by ([26]) oj < 1, so (1 — w) p ~ 3 is fine there. 
The same type of observation works for other cases above. We see that in all cases 
sign[(p — 1)$' — o;<&"] = — sign(p — 2). Therefore in the first two cases we have 
Z?2 < and this solution satisfies neither requirement (|13p nor requirement (|14|) . 
In the second two cases we have D\ > 0, and the function M can be a candidate 
either for M max or for M m i n depending on the sign of the second derivative My 3 y 3 . 
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Recall that (see (HDD ) 

In the case 32) we have sign[wi?2 + (p — l)R\] = — sign(p — 2), and therefore we need 
only to know signi?i = sign^' = sign 1) Since this solution is considered 

only in the sector p -r 2 -2/i < 2/2 < 2/1 (see Fig. [1]), we have 

G(yi + 2/2,2/1- 2/2) = (2yi) p - 1 \py 2 - (p - 2) yi ] > , (43) 

and id, being the unique positive solution of the equation 

G(u, 1) = (w + l) p_1 [w - p + 1] = -G( yi +2/2,2/1- 2/2) , (44) 

2/3 

satisfies the condition oj > p— 1. Therefore, signify = sign -^G{oj, 1) = signp(w + 
l) p ~ 2 (u; — p + 2) > 0, and so signiWy j, = — sign(p — 2), i.e., for p > 2 this is 
candidate for M max and for p < 2 this is candidate for M min . 

We are still considering the case 32). Recall that this function is defined not in 
the whole domain H + , but only in the sector ^-^yi < 2/2 < 2/1- To get a solution 
everywhere we need to "glue" this solution with that we obtained considering the 
case 2) (see (fTHD): 

M(y) = (2/1 + 2/ 2 ) p + C(y 3 - (2/1 - y 2 ) p ) • (45) 

To glue this solution along the plane 2/2 = ^rVi with that we just obtained, let 
us require from the resulting function to be continuous everywhere. From (|44p 
and (|43p we see that G(w, 1) = on this plane. Therefore, u = p — 1 and M = 
u; p t/3 = (p — l) p 2/3- The same value has solution (|45p on this plane for C = (p— l) p . 

Now we need to check that we get correct continuation in the sense that if the 
solution satisfies (|13p . then its continuation satisfies the same condition as well, if 
the solution satisfies (|14p . then the same is true for its continuation. The Hessian 
determinants will have the right sign automatically (actually D2 = identically). 
We need only to check the sign of 

M ViVi = My 2V2 = p(p - l)(( yi + y 2 r 2 - (p - i) p (yi - y 2 ) p - 2 ) 

in the domain — y\ < 2/2 < P_ ^r2/i, or m the initial coordinates < X2 < (p — 1)^1- 
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For p > 2 we have 

(z/i + y 2 r 2 = *r 2 < (p - ir m~ 2 

<(p-i)%r 2 = ( P -if( yi -y 2 r 2 , 

and for p < 2 we have 

(yi + 2 = *r 2 > (p - i)^- 2 xr 2 

>(p-i)%r 2 = ( P -if( yi -y 2 r 2 . 

This means that M is a candidate for M max if p > 2 and a candidate for M m ; n if 
p < 2, as it has to be. 



Let us rewrite expression (|45p in the same form, as it was made in 

M-Cy 3 = ( yi + y 2 f - C(yi - y 2 ) p = x\ - Cx\ . (46) 

Therefore, if we change a bit the definition of G defining it on the quadrant z% > 
as follows 



G p (zi,z 2 ) 



if zi<(p- l)z 2 , 



_ ( Zl + za)^ 1 [z! - (p - l)z 2 ] , if zi > (p - l)z 2 , 
then we can write two our solutions M on H + in an implicit form as before: 

G{yi + y 2 , y\ - y 2 ) = y-iG(uj, 1) . 

or solutions B on f2 + 

G(x2,ari) = G(S*,x|), 



(47) 



(48) 



In the case A 2 ) we again consider exactly the same G p from (|47p . In a similar 
way we can glue continuously the solution in case A 2 ) found in the sector —y\ < 



V2 < ^fyi 



G(l,u) = (u + - (p - = -G( Vl - y 2 , Vl + y 2 ) , (49) 

V3 



which is the same as 



G( Xl ,x 2 ) = G(xg,B p ] 



(50) 
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with the solution (|45p along the line 2/2 = Here we have to take C = (p' — l) p , 

because on the line 7/2 = ^-jpVl we have G(1,lj) = 0, i.e., u = p' — 1. Now, in the 
sector — 2/1 < 2/2 < ^jp2/i we have 

pu p - 2 R?H 2 p-2 



1/32/3 2/3 w + 1 ' 

Therefore, signM^ w = sign(p — 2), i.e., for p < 2 this is candidate for M max and 
for p > 2 this is candidate for M m i n . 

In the "dual" sector X2 > (p' — l)xi (or 7/2 > 2 ~^2/i) for p > 2 we have 



(2/1 + y 2 r 2 = *r 2 > &/ - iy a *r a 

> (p>-l)P X P- 2 = (p>-lf(y 1 -y 2 y-\ 

and for p < 2 we have 

(2/1 + y 2 r 2 = *r 2 < & - i) p - 2 ^r 2 

< ( P ' - i)%r 2 = o>' - i)f( yi - 2/ 2 ) p - 2 . 

This means that M is a candidate for M max if p < 2 and a candidate for M m j n if 
p > 2. 

Using the same "generalized" definition (|47p of the function G we can write our 
solutions M on S + in an implicit form as before: 

G{yi - 2/2, 2/1 + 2/2) = V3G(1, w) • 

or solutions -B on f2 + 

G(xi,x 2 ) = G(x|,S ? ) J (51) 

which should give, as we said above, the candidate for B max for p < 2 and B m j n for 
p > 2. Notice that for p > 2 the candidate for, say, B max is given by equation (|4"8]h 

It is a bit inconvenient to use one equation for, say, B max if p > 2 (this will 
be (j48j) ). and another one (this will be ([ST]) ) for the same B max if p < 2. We note 
that after interchanging role of Zi and replacing p by p' we get the scalar multiple of 
the original expression in both lines of (I47h . This allows us to give one expression 
for i? max for all p using notation of p* = max{p,p'}. In such a way we come to 
formula for F p , where we introduce additional scalar coefficients to make this 
function not only continuous but C 1 -smooth everywhere in f2 + . This smoothness 
guarantee us that the solution B is C 1 -smooth as well. 
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5. Proof of Theorem [TJ Verification theorem. 



Exactly in the spirit of Stochastic Optimal Control theory we wrote the PDE (|15p , 
we solved it in the previous section by building B which solves the equations of 
Theorem [1] (these are the same equations as (|48p , (|50p ) . Now continuing in the 
spirit of general results of Stochastic Optimal Control theory [FR] , |WF| we need to 
prove that these solutions in fact are equal to B max , B m i n . In Stochastic Optimal 
Control theory such proofs are called verification theorems, and they state roughly 
that if the solutions have a certain smoothness (often even slightly less than C 2 ), 
and if the domain is convex, then we are fine. 

From now on we denote by -B max the unique positive solution of the equation 

F(|x2|, |xi|) = F(B P , £3 ) and by B m \ n the unique positive solution of the equation 
I I 

F(|xi|, \x2\) = F(x^ , B p ), where the function F = F p is defined in fl3J). Existence 
and uniqueness of the solution follows from the fact that F(zi,z 2 ) is strictly in- 
creasing in Z\ from — — Xfz\ till +oo as Z\ runs from to +oo and it 
is strictly decreasing in Z2 from p(p* — Vf~ x z^ till — oo as Z2 runs from to +oo. 
Indeed, the first partial derivatives of F are 



p-i 



p(l 



(zi + z 2 y- 2 [pz l 



if z\ < (p* - l)z 2 , 

{(p-l)(p* -l)-l)z 2 ] 
if zi > (p* - l)z 2 ; 



(52) 



-{p* - Ifpzl 



if zx < {p* - l)z 2 



-p{l 



_L\p-l/ 



(53) 



[ ( Zl + z 2 y~ 2 [( P * - P ) Zl + P (p* - i) Z2 ] , 

if zi > (p* - l)z 2 . 

Note that both derivatives are continuous everywhere (even at the origin, where 
they vanish). Moreover, F Zl > if z\ > and F Z2 < if z 2 > 0, i.e., F is strictly 
increasing in z\ and strictly decreasing in z 2 . 

In the case of S max we look for a solution of the equation 

F(B*,xl) = F(\x 2 \,\xi\) 



or 



F(w,l) = —F(\x 2 \,\xi\ 

X3 



BURKHOLDER'S FUNCTION 



27 



Thus, we get a continuous solution u(x) everywhere except the plane X3 = 0, 

where w is not defined. But we can easily estimate the behavior of w nearly the 

1. 

line X3 = x\ = 0. Since F is decreasing in z 2 and < |xi| < x% , we have 



— , 1) < F(uj, 1) = F{-^, — r^) < F{—^^ 

3 



•j- o 3 3 

Since F is increasing in z\, we get 

\x 2 \ 

—if P < w < w o , 
where wo is the solution of the equation 

^o + lf-Vo-^ + lH — 

^3 

Whence ojq > p* — I and 



(wo - p* + l) p < (w + l)^ 1 (wo - V* + 1) 



F2 



^3 

i.e., 

Wo < p - 1 + 



Mp ' 

■"3 

Therefore, for B = w p £3 we have the following estimate 

\x 2 \ p <b< (is2i+(p*-i)4 / y, 

which gives the continuity near X3 = 0. Thus, the solution -B max is continuous in 
the closed domain fi. 

Similar considerations gives us the continuity of B mm . In that case we have the 
equation 

F(l,w) = — F(\ Xl \,\x 2 \) , 

X3 



and hence 



X^ X 

Now F(l,w) is decreasing in w, therefore, 

\x 2 



< W < Wq , 



where wo is the solution of the equation 

1 - (p* - i)M = -(p* - ir 



^3 
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i.e., Wq = (p* — 1) v + |^2| p /x3 and for B = u p xs we have the following estimate 

\x 2 \ p <B< \x 2 \p + (p* -l)~ p x 3 , 

which gives the continuity near X3 = 0. Thus, the solution B m { n is continuous in 
the closed domain as well. 

First step of the proof is to check that the the main inequality (concavity (jlip 
for the candidate -B max and convexity f)12|) for the candidate -B m in) is fulfilled if the 
points x + ,x~ satisfy the extra condition on their coordinates: 

\ x t ~ X l\ = \ X 2 ~ X 2 I • (^) 

This was almost done in the preceding section, when constructing these candidates. 
We know that the Hessians of our candidates have the required signs everywhere 
in our convex domain f2 except, possibly, the planes x\ = 0, x 2 = 0, and, either 
\ x 2\ = {p* — l)|^i| for -B max or |xi| = (p* — 1) | a^2 1 for B m \ a . On these hyperplanes our 
solutions are not C 2 -smooth, but this does not prevent them from being correctly 
convex (for the 32), p > 2 and 4 2 ), p < 2 cases) and correctly concave for the rest 
of the cases (namely, for the 3 2 ), p < 2 and 4 2 ), p > 2 cases). This one checks just 
by calculating directly the sign of the jump of the derivative. Namely, one fixes the 
line Lt = a + bt in the direction of the vector b = (bi,b 2 ,b 3 ) such that |&i| = \b 2 \. 
We need to prove the concavity of B, the candidate for B max , and the convexity of 
B, the candidate for B m ; n on Lt. At any point of Lt, which is not the intersection 
of Lt with the aforementioned hyperplanes, this concavity (convexity) follows from 
the previous section, this is how the candidates for B max , B m j n were built in (|48p . 
(|50p . At the points of intersections of Lt with the hyperplanes one can check the 
sign of the jump of the derivative of B(a + tb). We leave this as an exercise for the 
reader. 

Let the triple of points satisfies the following relations 

I rf* I /y» I I fy I /y» I ry* I I ry* I IP ry* I rv* IP 

I 1 1 I — I 2 2 I ' 3 — — I 1 I ' 3 — — I 1 I ' 

--++-+ ± ( 55 ) 

x = a x + a x a + a = 1 , a > 0. 
Now we have our solution -B max the following main inequality (biconcavity) 

B{x) -oTB{x-) -a + B{x + ) > 0, (56) 
and the opposite main inequality (biconvexity) 

B{x)-a~B(x-)-a + B(x + )<0 (57) 
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is true for the solution B m \ n . 

Lemma 2. If a continuous in 0, function B satisfies the main inequality (|56p 
and the boundary restriction B(x\,X2, \ x\\ p ) > | X2\ p , then B > B max . If it satis- 
fies (j57|) and B(xi,X2, \ x\\ p ) < \ X2\ p , then B < B m j n . 

Proof. Let / = [0, 1] and J denote an arbitrary its dyadic subinterval. As always 
J+, J_ are two sons of J. Let us fix two bounded measurable test functions /, g on 
/ such that \(g,hj)\ = |(/, hj)\ for any J. Put 

The fact that \(g,hj)\ = |(/, hj)\ exactly guarantees that x + , x~ satisfy the as- 
sumptions of ([55]) and we can rewrite inequalities (|56|) and (|57|) with x = Xj in the 
form 

\J\B(x) -cT\J-\B(x-) - \J+\a + B(x + ) > 0, (58) 
\J\B(x) -or\J-\B{x~) - \J + \a + B(x + ) < 0. (59) 

Let J n denotes the set of dyadic subintervals of n-th generation, i.e., Jq = 
{/}, and J n is the set of suns of elements from J n -\. So, adding up all our 
inequalities (|56l) with x = Xj for J G J n -\ we get 

\AB{x 3 ) > Y, \J\B( Xj ). 

Adding up these inequality over n from 1 to iV we get 

B{x)> \J\B{xj) = f B(x N {t))dt, 
JeJ N Jo 

where x (t) is the step function equal to Xj if t E J for every J from Jn ■ 

Notice that {</?), —> <p(t) almost everywhere when runs over a family of nested 
intervals shrinking to the point t for an arbitrary summable function ip. Therefore, 
almost everywhere 

x N (t) = ((f) r , (g) r , ( \f\ p )j) -+ (f(t),g(t), \f(t)\ p ) as N -+ oo 

and since B is continuous, we have 

B(x N (t))^B(f(t),g(t),\f(t)\ p )>\g(t)\ p . 

Now using Lebesgue dominant convergence theorem we come to the estimate 
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for every pair of bounded measurable functions f,g. And finally approximating 
arbitrary f,g G L P {I) by its cut-off functions and using monotone convergence 
theorem we can extend this inequality to the set of arbitrary possible test functions 
/ and g what niGcins exactly thcit B max (x) < B(x). 

For the case of B rmn in all these considerations we need to change the sign of 
inequalities only, and we will get B mm (x) > B(x) for B satisfying (f57|) . 

□ 

We are left to prove the opposite inequalities 

B max ^ B max (x) and B m ; n < B m 'm(%)- 

This can be done by reversing the reasoning in the lemma above. Using the fact that 
domain £1 = {x = (xi,X2,X3) : X3 > |xi| p } is foliated by the straight line segments 
(extremal trajectories) it is possible to construct the sequence of test functions 
fn,9n corresponding any given point x G 0, and such that {\g n \ p )j — > B{x). This 
just supply us with the required inequality. The reader can see how this type of 
reasoning is done in [VaVo2| . The main idea is to travel along the extremal trajec- 
tories starting from x £ to build a net N := {x + , x~, x ++ , x H ,x , x 
All points of the net should belong to f2, and we put them on the same extremal 
trajectory on which x lies for a while. If one of them, say, z hits the boundary: 
d£l (parabola) we stop building children z + ,z~. But then one of them, say, ( can 
hit the special hyperplanes x\ = or X2 = 0. In this case we choose C + >C~ hi 
such a way that they lie in different quadrants very close to Then we start anew 
a building of the net for £ + and C~ separately. The closer £ + >C - are to C the 
smaller will be difference {\g n \ p ) I — B(x). In such a way for arbitrary e we obtain 
the inequalities 

B max (x) > -B max (x) - e and B min (x) < B min (x) + e . 

The reader can address to |VaVo2| to understand how the net M generates a 
required pair of functions /, g. But in the proof of the following lemma we only 
state the result of the described construction that supplies us with a recursive 
definition of / and g. 

Lemma 3. The functions -B max and B m { n satisfies the inequalities 
B m ax(x) > -B max (x) and B min (x) < -B min (x) . 
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Proof. We construct an extremal sequence of pairs /, g for the function B max (x) for 
p > 2 and for some point x on the plane x\ = 0. For p < 2 the same construction 
works for x<i = 0. For all other points we "glue" the extremal pairs form the known 
functions on the ends of the extremal trajectories. The detailed explanation how 
to do this can be found e.g. in |SVj . 

Take an arbitrarily small e and recursively define the following pair of test func- 
tions 



/(*) 



and 



g(t) 



-c, 
7/(£§) 
c, 

d-, 



< t < e, 
e < t < 1 
e < t < 1 



1 



< t < e, 
e < t < 1 - e , 

£ < t < 1. 



where the constants c, d±, and 7 will be defined from the conditions that guarantee 
that this is an admissible pair of test functions corresponding to a given point 
x = (0, X2,x%). It is not difficult to see that these formulas correctly define / and 
g almost everywhere on [0, 1]. 

It is evident that x\ = (/) = 0. Since 

(0>[o,i] = £ ( d ~ + d +) + ( X " ^Hci] > 
the condition (g) [0 x] = X2 is 

e(cL. + d+) = (1 - 7 + 2e 7 )x 2 . 

The condition (|/| p ) [0 x] = X3 gives the relation 

2ec P = (1 _ 7 p + 2e7 P )x 3 . 



(60) 

(61) 
1-e 

c+ __ = 1 __ f _^_ d+ , (62) 

The left point x + is already on the boundary <9f2 (the functions are constants on 
/ + ) and we have nothing to split. The left interval I~ we naturally split at the 



Two more relation we obtain by using condition \xf — x 1 
be the first splitting point then the condition xf — xj~ 

ec ed- + (1 — 26)7x2 



0C 9 



— x 2 I • Let t 

xt is 
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point t = e. Then the condition x± — x~[ = x\ — x 2 

c = 7x2 — d- . (63) 

From ([63]) and (j62|) we get 

d- = 7x2 - c , 

(64) 

d+ = 1X2 -l^2e- 
Now we can plug in (|64p into (|60p and obtain in result 

7 = 1 + 2e-^^- • — . (65) 
' l-2e x 2 

Let Co be the limit value of c as e tends to 0. (By the way, cq has clear geometrical 

meaning: this is the first coordinate of the end point on d£l of the extremal line 

the second end of which is our initial x.) Then 

7-1 + 2^, 

Y ~ 1 + 2 p^£ , 
X2 

l-(l-2e) 7 p ~2(l-^)e, 



X2 



and ([6T]) turns into equation for cq: 



4 = (1 - ^°)* 3 , (66) 

X2 

which evidently has unique solution cq, < cq < Further we get 

d± — > x 2 - c (67) 
and therefore we are able to write down in term of Co the average we interested in: 

_ (<£ + <)£ 2(X 2 - C )P _ X 3 (X 2 ~ Cp)P 

W\ hoAl l _ {l _2 £ ) 1 P ^2(1-^) eg • 

If we introduce 

_ 1 1 

then from (|68|) we get Co = • We can plug this expression in (fBHj) and conclude 
that cj is the unique solution of the equation 

x p 2 = (co + l) p ' 1 (u + l-p)x 3 , 

but this is just the equation F p (x2,0) = F p (uj, 1)3:3, whose solution by definition is 

'£ max (0,22,X 3 )\ I 



X3 
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Therefore, 



B mSkX (0,x 2 ,x 3 ) = lim( |ff| p ) [0il] , 



£->0 

what yields the desired inequality 

-Bmax(0, X 2 ,X 3 ) < B(0, X 2 , X 3 ) . 

How to obtain this inequality for arbitrary x £ SI was explained in the beginning of 
the proof. In the same paper [SV| the reader can find the detail explanation how 
to pass from one splitting to another one (e.g. to a dyadic family of intervals). 

□ 

6. Function u p from function B 

We found Burkholder's functions B max and B m j n as claimed in Theorem [TJ As 
a corollary we immediately we get the sharp constant in Burkholder's inequality: 

Theorem 4. Let I = [0, 1], (f) r = x\, (g) 7 = x 2 , g is a Martingale transform of 
f , and \x 2 \ <\x\\. Then 

< |#>, < (p* - l)W>i ■ 
The constant p* — 1, where p* := max(p, -^j) is sharp. 

Proof. We just analyze the form of function B max from Theorem[T]and immediately 
see that 

sup = yp — ly . 

x&Cl,\xi\<\xx\ ^3 

□ 

Theorem 5. Let L = [0, 1], (f) I = x\, (g)j = x 2 , g is a Martingale transform of 
f, and \x 2 \ < \xi\. Then 

(im,<(p*-i) p (i5a. 

The constant p* — 1, where p* := max(p, -^-) is sharp. 

Proof. We just analyze the form of function B m ; n from Theorem [T] and immediately 
see that 

B m in (x 1 ,x 2 ,x 3 ) 



inf " r """ iJ = (p*-l 

x&Q, |a!2|>|xi| 2:3 



□ 
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Remark. The same analysis shows that ( \g\ p )j < (p* — l) p { \f\ p )j if and only if 
I x 2 \ < (p* - 1)| xi\ in Theorem HI and in TheoremE] ( < (p* - l) p ( |5| p ) 7 if 
and only if \x 2 \ > (p* — 

Notation. Below we use /3 p := (p* - Put 

<Amax0i,£2) := sup [B max (a;i , x 2 , x 3 ) - (3 p x 3 ] , 
4>mm(x 1 ,x 2 ) := inf [B min (xi, x 2 , x 3 ) - P~ l x 3 ] . 

x 3 :(xi,X2,X3)£Q 

These functions are defined on the whole R 2 . 
Definition. If for all pairs of points x^ 6 R 2 such that 

\x~± — x± \ = \x% — % 2 \ an< ^ x = ~n (^9) 

the function <j) on R 2 satisfies the condition 

<P(x)- HX ~ )+ 2 * iX+) >0, (70) 
then it is called zigzag concave. If the opposite inequality holds 

«») - ^l±^j < (71) 

the function <fi is called zigzag convex. The next theorem gives an independent 
description of (p max and 4>min- 

Theorem 6. Function max is the least zigzag concave majorant of the function 
/i m ax(^) := | x 2 \ p — (3 p \ xi\ p . Function (p m m is the greatest zigzag convex minorant 
of the function h m - m (x) := | x 2 \ p — fi p 1 \ x\\ p . 

Remark. Notice that this is slightly counterintuitive: B max (x) — f3 p x 3 is zigzag 
concave for any fixed x 3 , and the supremum of concave functions is not usually 
concave. The same is true about infimum of convex functions. 

Proof. Let x^ and x are as in [69j It is obvious that 4> max is zigzag concave. One 
verifies this just by definition. In fact, if for any x~ G R 2 we can choose such that 
the supreniuni in the definition of ^ max is almost attained, i.e., B max (x ) — f3 p x 3 > 
4>ma,x( x -) — e for a given e and the same for x + e R 2 . We define x 3 = X3 1^ 3 and 
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x = (x,x 3 ). Then using (|56|) we can write 
<f>max(x) > B max (x) - (3 p x 3 

. B max (x ) + B max (x^~) Xg -f- Xg 

- 2 ^ p 2 

0max(ff~) + ^max(^ + ) _ 

2 £ ' 
what yields (j70j) . Inequality (jTTj) is totally similar. 

As sup is bigger than lim we conclude 

<Pma, x ( x ) > lim [B max (x) - (3 p X 3 ] = | X 2 | P - P p \ Xi\ P = h max (x) . 

As inf is smaller than lim we get analogously 

</>min(z) < lim [B min (x) - P^xz] = I x 2 \ p - /3" 1 ! xi\ p = h min (x) . 
This is because the boundary values of B max and B min are | x 2 \ p - 

We are left to see that (/> max is the least such majorant (and a symmetric claim 
for m i n ). Let ip be a zigzag concave function such that 

0max > V> > ^max • (72) 

Consider function ^(x) := ip(x) + /3 p X3. It is immediate that \& satisfies ([56]) . On 
the boundary of Q we have ^(x) > | X2\ p , this is just by the right hand side of ([72]) . 
Then Lemma [2] yields 

> B max (x) . 

Then, obviously, 

^(x) = sup - /3 p x 3 ] > sup [B max (x) - [3 p x 3 ] = max (x) . 

So we proved that 4> m£LX is the least zigzag concave majorant of /i max . Symmetric 
consideration will bring us the fact that 4> m i n is the largest zigzag convex minorant 
of h m [ n . 

□ 

The reader should look now at function F p from Theorem [TJ It would be inter- 
esting to obtain the formulae for </> max and </> mm , especially using this F p . It would 
be also interesting to understand the role of function 

u p (x 1 ,x 2 ) := p{\ - ^) p - X (\x x \ + \x 2 \f-\\x 2 \ - (p* - 1)\ X1 \) , (73) 
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mentioned in the introduction and used repeatedly by Burkholder. May be it is 
equal to <^ max ? The answer is "no" , but we can prove the following 

Theorem 7. 

<f>max(xi,X2) = F p (\ X 2 \, | Xi\) . (74) 

Proof. We shall consider only the case p > 2, the case p < 2 is similar. Due to the 
symmetry with respect to the change of x\ to —x\ and x 2 to — x 2 it is enough to 
check equality (|74p in the quadrant x\ > 0, x 2 > 0. If x 2 < (p — l)x±, we get an 
explicit formula for B max from Theorem[TJ B max (x) = x\ + (p — l) p (x3 — x p ), and 
therefore, 

<Pmax{xi,x 2 ) = sup [B max (x) - (p- l) p x 3 ] =x p 2 -{p- l) p x\ = F 2 (x 2l xi) . 

So in the rest of the proof we shall consider only the domain {x = (x\,x 2 ) ■ < 
(p — l)x\ < x 2 }. Moreover, since both functions ^ m ax 

and F p are p- homogeneous, 

it is sufficient to check (l74j) on the interval S := {x: < px\ < l,x\ + x 2 = 1}. 
(Indeed, the condition px\ < 1 on the line x\ + x 2 = 1 means x 2 > (p — l)x\.) 

The function F p is linear on S: F p (l — x±, x\) = ^ P ~l-2 — (1 — px\). Now we check 
that 4> ma _ x is linear as well. To this end we check inequality 

, ( , 0max(a;i - a, x 2 + o) + (?Wx(xi + a,x 2 - a) 

( Pmax\X) <. ^ V / 

for all x £ S and sufficiently small a, this just means linearity of ^ max on 5, because 
the opposite inequality follows from the zigzag concavity of </>max 

Fix x £ S and e > 0. Take X3 such that S(x) — /3 p X3 > 4> mSLX (x) — e. Due to 
condition x 2 > (p — l)xi the extremal trajectory L x of B max passing through the 
point x = (x, X3) is not vertical, it hits at some point the plane x\ = 0. Therefore 
we can take two different points x^ = (x^,x^) on L x such that x = \{x + + x~). 
We know three things: 

B max (x) - f3 p x 3 > 4> ma , x (x) - e , 

B m ax (x ) ftpX^ ^ </ ) max(2' ) j 
B ma x(^ ) fipXr^ ^ ^max^ ) • 
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Since the function B max is linear along L x , we can write the following chain of 
inequalities 

(x) - e < B max (x) - f3 p x 3 

[B max (x+) - PpXz] + [B max (x") - PpX^] 
~ 2 

(s+) + (j) max (x-) 

2 

Since e is arbitrary, we come to the desired convexity (|75p . 

Function F p (x2, x\) is a concave C 1 -smooth function majorazing /i max on S\ This 
is immediate from its formula. Functions </>max (xijX?) and F p (x2,xi) are linear on 
S and at the point x = x p =: ( | , 1 — ^) both are equal /i m ax(%>) = 0. Therefore, 
to prove that they are identical it is sufficient to check that their derivatives at x p 
along S are equal as well. Since </> ma x is a majorant of /i max and both functions are 
equal at x p , then the left derivative of </> max at x p is not grater than the derivative 
of h max at this point. On the other hand, since 4> max is the least majorant it is 
not grater then F p , i.e., its left derivative at x p is not less that the derivative of F p 
there, but latter coincides with the derivative of h max . Hence all three derivatives 
along S are equal at the point x p and we proved (f> m3X (xi, X2) = F p {x2,xi). 

□ 

Theorem 8. 

<th^{xx,x 2 ) = -ip* - iy p F p (\ xi\,\ s 2 |) • (76) 
The proof of this Theorem is absolutely similar to the proof of Theorem [71 

Burkholder often used function u p from (|73p. To demystify it let us notice 
that it is also p-homogeneous and as such can be considered only on the segment 
x\ + X2 = 1, Xi > 0. On this segment function u p becomes linear. It is a majorant 
of hp, and its graph is tangent to the graph h p exactly at point x p on S, where h p 
vanishes. It is not the least zigzag concave function greater than h p (of course not, 
0max is such), but it is the least zigzag concave function larger than h p and such 
that on all segments {x: Xi > 0, x\ + X2 = const} it is not only concave, but also 
linear. (Keeping in mind the symmetries x — > —xi,X2 — > —X2 we can consider the 
first quadrant only.) 

This is already proved, and we leave the detailed reasoning to the reader. One 
more thing we want to mention is that we could have considered a slightly more 
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general problem. Namely, instead of majorazing the function /i max (xi, x 2 ) = \ x 2 \ p — 
(p* — l) p \ xi\ p we could have started with any function 

h c (xi,x 2 ) ■= \x 2 \ p - c\xi\ p . 

The reader can easily see that we have proved the following theorem (of course 
Burkholder already proved the most of it long ago). 

Corollary 9. The smallest c for which there exists a zigzag concave function <p c ma- 
jorazing h c is equal to (p* — l) p . For this c the least zigzag majorant is F p (\ x\\, \ x 2 \). 
The smallest c or which there exists a zigzag concave function <p c majorazing h c such 
that it is linear on the segment {x: Xi > 0, x\ + x 2 = const}, symmetric and p- 
homogeneous is equal to (p* — l) p . For this c the least zigzag majorant linear on 
{x: Xi > 0, xi + x 2 = const}, symmetric and p -homogeneous is u p (xi,x 2 ). 

Remark. Notice an interesting thing which we do not know how to explain. Given 
function B max from Theorem Q] we can easily diminish the number of variables and 
construct (/> max . But amazingly we can also find B max if only ^ max is given. In fact, 
Theorem [7] gives the formula for max via F p . Then F p allows us to find B max . If 
we now combine Theorem [T] and Theorems [7HS1 to conclude that 

Corollary 10. Given a point x & SI, if we know ^ max we can find B max (x) by 
solving equation: 

1 1 

0max(3'l > X 2 ) = ^max(2'3 ) Bmax) • 

Symmetric formula allows to find B m i n if 4> m i n is known: 

1 1 

(t>min{x 2 ,X 1 ) = min (B^ in ,x|) . 
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